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Abstract We consider a one-dimensional gas of positive and negative unit 
charges interacting via a logarithmic potential, which is in thermal equi- 
librium at the (dimensionless) inverse temperature j3. In a previous paper 
[Samaj, L.: J. Stat. Phys. 105, 173-191 (2001)], the exact thermodynamics 
of the unrestricted log-gas of pointlike charges was obtained using an equiv- 
alence with a (1 + l)-dimensional boundary sine-Gordon model. The present 
aim is to extend the exact study of the thermodynamics to the log-gas on a 
line with alternating ± charges. The formula for the ordered grand partition 
function is obtained by using the exact results of the Thermodynamic Bethe 
ansatz. The complete thermodynamics of the ordered log-gas with pointlike 
charges is checked by a small-/? expansion and at the collapse point /3 C = 1. 
The inclusion of a small hard core around particles permits us to go beyond 
the collapse point. The differences between the unconstrained and ordered 
versions of the log-gas are pointed out. 

Keywords Two-component log-gas • Charge ordering • Exact thermody- 
namics • Thermodynamic Bethe ansatz 



1 Introduction 

We study thermal equilibrium properties of a symmetric two-component 
plasma (Coulomb gas) which consists of mobile pointlike positive and nega- 
tive unit charges ±1, confined in a two-dimensional (2D) domain of points 
r = {x,y). The charges are immersed in a homogeneous medium of dielec- 
tric constant 1, the system as a whole is electroneutral. In Gauss units, the 
Coulomb potential </>(r), induced by a unit charge at the origin 0, is given by 
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the 2D Poisson equation 



27r5(r). 



(1.1) 



Such definition of the 2D Coulomb potential maintains many generic prop- 
erties of "realistic" three-dimensional Coulomb systems with 1/r potential, 
e.g. the screening sum rules |25j . In an infinite space, the solution of Eq. (II. 1|) 
reads 



where r = |r| and r$ is a length scale which fixes the zero point of the poten- 
tial. The interaction energy of two charges q and q' at the respective spatial 
positions r and r' is equal to v qq >(r, r') = qq'<j)(\r — r'|). The corresponding 
Boltzmann factor at the (dimensionless) inverse temperature /3 = l/{k^,T) 
reads exp[— /3v qq i(r, r')] = |r — r'|^ OT . For two oppositely charged pointlike 
species at distance r, the Boltzmann factor is integrable at small dis- 
tances r — > for small enough ft, ft < f3 c . Thus, the thermodynamics of 
pointlikc charges on a continuum space is well defined only in the high- 
temperature region ft < f3 c . For low temperatures /3 > /3 C , the short-distance 
(ultraviolet) collapse of oppositely charged point particles makes the thermo- 
dynamics unstable; the thermodynamic stability is restored by considering 
a short-distance regularization of the Coulomb potential, i.e. via hard cores 
attached to particles which prevent them from touching one another. The 
value of the collapse (inverse) temperature /3 C depends on the dimensionality 
of the compact domain A to which the charge system is constrained. 

In the case of a 2D domain A, the 2D integral J A d 2 r \r\~@ is finite at 
small distances provided that /3 < 2, i.e. f3 c = 2; there is no problem at large 
\r\ because the interaction is screened by the conducting system. A typical 
domain of this kind is an infinite space, A — > M 2 , which defines the bulk 2D 
Coulomb gas. The complete thermodynamics of the system of pointlike ± 
charges in the stability region j3 < 2 was derived in [33] ; for a short review of 
exact results for bulk and surface thermodynamics, together with asymptotic 
large-distance behavior of charge and density correlation functions, see [28] . 
The derivation was based on an equivalence between the 2D Coulomb gas and 
the (l+l)-dimensional sine-Gordon theory, transferring from that integrable 
theory some exact results obtained within the Thermodynamic Bethe ansatz 
(TBA). The extension of the exact thermodynamics of the bulk Coulomb 
gas beyond the collapse point f3 c = 2, with a hard-core regularization of the 
Coulomb interaction and in the region of low particle densities, was proposed 
in |23] . An electroneutrality sum rule and the leading short-distance behavior 
of pair correlation functions was used to go up to f3 — 3. Applying a system- 
atic short-distance expansion of correlation functions [21] , Tellez showed [3"2"] 
that one can proceed, in principle, up to the Kosterlitz-Thouless transition of 
an infinite order from a high-temperature conductor phase (a non-zero frac- 
tion of the positive and negative charges is dissociated) to a low-temperature 
insulator phase (the positive and negative charges form neutral pairs) which 
takes place at /?kt = 4 for low densities. 

In the case of a one-dimensional (ID) domain A, say the infinite line 
x G (— oo,oo), the ID integral J A dxx~ /3 is finite at small distances provided 




(1.2) 
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that j3 < 1, i.e. /3 C = 1. Such systems, to which we refer as log-gases, have 
evoked much of interest because of their relationship to various models of 
condensed matter. There exist two basic versions of two-component log-gases: 
without and with space restriction on the ordering of ± charges. 

— Unconstrained log-gas: The "standard" ID log-gas without any space 
restriction on the ± charges is related to dissipative quantum mechan- 
ics [29 , 5 , 20 , 6 and to the problem of non-equilibrium quantum transport 
through a point contact in a ID Luttinger liquid 24,11 . The lattice ver- 
sion of the model, which represents a kind of short-distance regularization 
of the Coulomb potential, was exactly solved (the grand partition function 
and the particle correlation functions) at (3 = 1 [TS] and /3 = 2,4 [HAfT?] . 
The conductor- insulator phase diagram was conjectured in [18]; various 
approaches indicate that the Kosterlitz-Thouless phase transition should 
occur at /3kt — 2, independently of the particle density. Thermodynamic 
properties of the continuous version of the model with pointlike charges 
was solved in the whole stability region /3 < 1 in [27], by exploring the 
TBA results for an equivalent (1 + l)-dimensional boundary sine-Gordon 
theory [TUIITTllII ig] ■ 

— Log-gas with charge ordering: The ID log-gas with ± charges required 
to alternate in space is equivalent to the Kondo problem with spin-^ 
impurity [2, 3, 30]. The lattice version of the system is exactly solvable at 
the "collapse" isotherm (3 = 1 [T3] . The charge ordered system is expected 
to exhibit a dielectric phase at an arbitrary temperature |14) . 

The asymptotic large-distance behavior of the particle correlation functions 
was studied in detail for both versions of the ID two-component log-gas in 

The aim of the present work is to extend the exact study of the thermo- 
dynamics to the two-component log-gas on a line with alternating ± charges. 
The basic formula for the ordered grand partition function is obtained by 
using the TBA results [T2"l[T3~] , associated with a fusion relation between the 
grand partition functions of the unconstrained and ordered log-gases. The 
complete thermodynamics of the ordered log-gas with pointlike charges is 
derived. The results are checked by a small-/? expansion and at the collapse 
point /3 C = 1. The inclusion of small hard cores around particles permits us to 
extend the thermal analysis beyond the collapse point j3 c = 1. The important 
differences between the unconstrained and ordered versions of the log-gas are 
pointed out. 

The paper is organized as follows. 

We start with a brief recapitulation and an extension of the results for 
the unconstrained ID log-gas [22] in Sect. 2. Sect. 2.1 deals with the case of 
pointlike charges. The inclusion of a small hard core around particles, which 
permits us to study the thermodynamics beyond the collapse point j3 c = 1, 
is the subject of Sect. 2.2. Here, we apply two methods: the one based on a 
perfect screening sum rule valid for conducting systems and the other, more 
general, based on the explicit definition of the grand partition function. 

The thermodynamics of the ordered ID log-gas is derived in Sect. 3. As 
in the previous unconstrained case, the system of pointlike charges is treated 
in Sect. 3.1 and the inclusion of the hard core to particles, via the definition 
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of the grand partition function, is worked out in Sect. 3.2. The differences 
between the unconstrained and ordered versions of the log-gas are pointed 
out. 

Sect. 4 is Conclusion. 

Auxiliary calculations are shifted aside to Appendices. A fusion relation 
between the grand partition functions of the unconstrained and ordered log- 
gases is used to rederive in an alternative way a relationship between the 
corresponding bulk pressures in Appendix A. The small-/3 expansion of the 
obtained density-fugacity relationship for the ordered log-gas is checked by 
microscopic calculations in Appendix B. 



2 Extension of the results for the unconstrained log-gas 

First we recapitulate and extend the exact results for the thermodynamics of 
the unconstrained log-gas obtained in [27]. We start with pointlike charges. 
The inclusion of hard cores to particles, which will allow us to pass through 
the collapse point /3 C = 1, is described in the subsequent part. 



2.1 Pointlike particles 

We go to an infinite line through the circle of radius R and circumference 
L = 2ttR, taking at the end the limit R — >■ oo. The position of a particle on 
the circle is specified by the angle ip £ [0,27r). Since the distance between 
particles at angle positions ip and ip' is 2R\ sm[(<p — ip')/2]\, the 2D Coulomb 
potential takes the form 



2R 



sin 



<P~<P 



(2.1) 



Defining the line coordinate as x = Rip, the ID realization of the potential 
(|1.2[) . 4>(x, x 1 ) = — ln(|x — a;'|/ro), results from this expression as the R — > oo 
limit. The Boltzmann factor of two charges q and q' at the respective angle 
positions ip and ip' reads as 



exp[-Pqq'<j>(<p,tp')] 



2R . ftp-tp' 
sin 



Pqq' 



(2.2) 



We shall work in the grand canonical ensemble. Due to the charge ± 
symmetry, chemical potentials of the species can be taken equivalent, /i+ = 



/i. Only neutral configurations with the particle numbers 



N- 



N are considered. The position angles of (+) charges will be denoted by tpi 
and those of (— ) charges by p[ (i = 1, . . . , N). The grand partition function 
is defined as 



N=0 



cxp(2^AQ 
(Nl) 2 
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ri(i<j)=i 
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2^ sin 









(2.3) 



where A is the de Broglie wavelength and the dependence on the inverse 
temperature (3 is omitted. We introduce the rescaled fugacity z + = z_ = z, 



J/* 



z = exp(/3/i)-2; 



A 



(2.4) 



which has the dimension [length] 2 1 . Factorizing out quantities of nonzero 
dimensions in (I2.3p . we end up with the representation 



S L (z) = 1 



E 

N=l 



(2tt)^ 2 zL 1 - 



where the dimensionless configuration integrals 



I2N — 



1 



dpi 
2ir 



(TV!) 2 

involve the Boltzmann weights 



2tt 



dtp 



N 



2n 



2 A 



'27V, 



(2.5) 



2tt 



S 2 jv({^},{^}) 



(2.6) 



B 2 at({^}, {¥?':}) = 



ri(I<j)=i 


2 sin ( 


2 /. 




2 sin 
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(2.7) 



In the large- L limit, we pass from a finite circle to an infinite ID line, x € 
(—00,00). According to elementary thermodynamics, the quantity \o.EEl{z) 
is expected to be extensive. In particular, the bulk pressure of the charge 
system P(z), defined by 



(3P(z) = lim -lnS L {z), 

L— >OQ Li 



(2- 



has a well-defined finite value. The only expansion parameter in the series rep- 
resentation (I2.5[) is the dimensionless combination zL x ~^ , so that (3P{z) oc 
z 2 -p . The proportionality /3-dependent factor was found in |27j by using the 
TBA results obtained for an equivalent boundary sine-Gordon model [TU1ITT1 
Hl[5] . We take the TBA results from [TO) ; the notation in that work is related 
to ours as follows: g — (3/2, t = 1/g = 2/(3, q — exp(ing) — exp(i7r/3/2), 
T = l/L and a scale Tb = The expansion parameter in (|2.5I) was 

expressed in terms of the length scale Lb as 



z(2nf/ 2 L 1 -^ = r 



L 



r 



2^r 





2(2-/3) 



)J 



(2.9) 
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where r is the Gamma function [19] , For large L, the grand partition function 
is given by 



L— >oo V 2 



L 



1 



L 



b 2< 



2(2-/3); 



(2.10) 



Eliminating Lb from the couple of equations and using formulas for the 
Gamma functions [19] 



r(i-y)r( y ) 



sin(7ry) 

with real y, we finally arrive at the bulk pressure 



r| rf(^t4) 



1 ,r(l=4\r 



27r 3/2 \2-(3 



2(2-/3) 



2?rz 



The one-body densities are defined as thermal averages 



(2.11) 



(2.12) 



n q (x) = n q 



\ j 



6 qig Mx-Xj) ) , 



(2.13) 



where 5 qtq i is the Kronecker symbol, S(x — x') the Dirac delta function and 
the index j runs over all particles. For the considered charge ± symmetry, 
the species densities are equal to one another: n + = n_ = n/2, where the 
total number density of particles n is given by 



n(z) = z-^-pP(z). 
oz 

Thus the density-fugacity relationship reads as 
1 „ (\-B\ „ / P 



n 



-ril=± )r 



2?r 



r(/3/2) 



The small-/? expansion of the rhs of this formula 

-4- = 2/3^ exp ([C + ln(27r)] £ + 0(/3 2 ; 

7^? V 2 



(2.14) 



(2.15) 



(2.16) 



where C is the Euler constant, was checked in Appendix of [27] by using a 
renormalized Mayer expansion. Note that this expansion is non-analytic due 

to the appearance of the term /3^F. On the other hand, the series in the 
exponential is analytic in /3. 

For a fixed z and in the limit (i -> l - , the term r((l — /3)/(2 — /3)) ~ 
1/(1-/3) implies that the particle density n exhibits the expected collapse 
singularity 

4z 2 

" ~ ( 2 - 17 ) 
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This singular behavior can be deduced indirectly by using a perfect screening 
sum rule for the one-body densities of the conducting system [25] . 



dx[U q ,- q (x)-U q>q (x)}, (2.18) 
where the Ursell functions are defined by 

U q<q >(x,x') = U q>q *{\x - x'\) = ( y^,$ q , q:i S(x - Xj)6 q/ . qk S(x' - x k ) \ - n q n ql . 

w* / 

. (2 - 19) 

The short-distance behavior of the Ursell function for a positive- negative pair 
of charges is given by the Boltzmann factor of the pair interaction [2"2"ll21| . 

U q ,- q (x) ~ z 2 \x\-P. (2.20) 

x— >0 

For /3 — > 1~, the integral in (|2 . 18[) is dominated by this short-distance be- 
havior and we have 

-2~f_^W = ^" = T^ + 0m - (2 ' 21) 

where £ is a screening length of the Coulomb potential. This derivation of 
the singular behavior (|2.17p points out the two-body nature of the collapse 
phenomenon. 



100 
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Fig. 1 The density-fugacity relationship versus the inverse temperature f3 for the 
unconstrained log-gas. The packing fraction rj = no, where a is the hard-core 
diameter. The solid curve describes pointlike charges (77 = 0). The dashed and 
dotted curves correspond to rj = 10~ 3 and 10~ 2 , respectively. 
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The plot of the density-fugacity relationship for pointlike particles (|2.15[) , 
well defined up to f3 c = 1, is represented in Fig. 1 by the solid curve. Although 
the function looks like monotonously increasing, this is not true: as is seen 
in Fig. 2, the function decreases and reaches the minimum for small f3. 
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P 

Fig. 2 Fig. 1 on a smaller scale to document the non-monotonous behavior of the 
density-fugacity plot for small /?. 

The virial equation of state, which relates the pressure P and the particle 
density n, takes a simple form 

0P = (l - f) n (2.22) 

which is characteristic for log-gas systems. 

Having the density-fugacity relationship, the full thermodynamics is ob- 
tained by passing from the grand-canonical to the canonical ensemble via the 
Legendre transformation 

F L (0,N) = ~hi3 L + nN, (2.23) 

where Fl is the Helmholtz free energy and N — nL. The (dimensionless) 
specific free energy, defined as / = (3Fl/N, thus reads 

f{p,n) = ln(An) - | ln(2r n) + \ (\ - Inn 
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The excess (i.e. over ideal) internal energy per particle and specific heat at 
constant volume are given by 



respectively. For the specific heat, we obtain explicitly 



fcr- 



P 2 



i>' 



, fl-P 



4/ 



P 



2(2 -Pf \J \2-P 
P 2 

2(2 - P) 4 

where the psi function and its derivative are given by 



2(2-/3) 



d °° 1 

4 {x ) = -lnr(x), ^(^g^-— 



(2.25) 



(2.26) 



(2.27) 



The specific heat exhibits the Laurent series expansion around the collapse 
point P — > 1 — , 

£ = ct-juV 0(1) - (2 - 28) 



2.2 Particles with hard cores 

To go beyond the collapse point P c = 1, we have to consider a short-distance 
regularization of the interaction Coulomb potential, especially for oppositely 
charged particles. We introduce a hard core of diameter a around each par- 
ticle which does not allow the charges to touch one another: 

{—qq'\n(\x — x'\/r Q ) if \x — x' | > er, 
oo if \x — x'\ < er. 

The total particle density now depends also on er, n(z,a). The hard-core 
diameter a, or more precisely the dimensionless combination 

£ = az^, (2.30) 

will be assumed to be small, £ — ?► 0, and we shall look for the leading hard-core 
correction to the pointlike system. 

We first derive the leading hard-core correction in analogy with the deriva- 
tion for the 2D Coulomb gas presented in [23] ■ The unconstrained log-gas 
system is in its conducting phase up to /?kt = 2, so that the screening sum 
rule (|2.18[) holds and the Ursell functions are well defined also beyond the 
collapse point. The difference U q ,- q — U q , q vanishes inside the hard core for 
the potential (|2.29j) . and the total particle number density is given by 



n(z, a) — 2 / dx [U q - q (x; z, a) - U q . q (x; z, a)} 

J a 



(2.31) 
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We can write 

Uq tg /(x;z,a) = U q>q >(x;z,0) + A M > (x; z, a), 



X > (J, 



(2.32) 



which defines A q>q / (x; z, a), vanishing when a — ¥ 0, as the change of the Ursell 
function due to the introduction of the hard core a to pointlike particles. 
Subtraction of Eq. (|2.3ip with a > and the same equation with a = leads 
to 



n(z, a) — n(z, 0) = — 2 / da; {U q .- q (x; z, 0) — U M (x; z, 0)] 



/•OO 

2 1 dx [A q - q (x; z, a) - A qiq (x; z, a 

J a 



)]• (2-33) 



We make a heuristic assumption: in the small-cr limit, we can neglect the 
quantities A q: ± q in (|2.33[) as they contribute to subleading hard-core correc- 
tions. Such assumption has already been made in [7|l8ll23] for the 2D Coulomb 
gas. Consequently, 



i(z, a) = n(z, 0) - 2 f dx [U q ,- q (x\ z, 0) - U q . q {x; z, 0)] . (2 
Jo 



34) 



As before, the integral is dominated by the short-distance behavior of the 
Ursell function (|2.20p and we arrive at the basic result 



n(z, a) — n(z, 0) — 4z z 



1-/?' 



(2.35) 



Strictly speaking, this result was derived in the stability region < (3 < 1 
where n(z,0) is well defined by (|2.15p . It is however reasonable to assume 
that the formula (|2.35[) can be analytically continued beyond /3 C = 1 because 
both the sum rule (|2.31|) and the leading short-distance behavior (|2.20l) . 
which were crucial in its derivation, remain valid up to /?kt = 2. 

To provide a convenient representation of (|2.35D , we express the dcnsity- 
fugacity relation for pointlike particles (I2.15|) in the form 



n (z,0) = t-^-ttZ— . 



where the function 
1 



4^-3/2 



r 



3-2/3 



r 



1-/3' 
13 



2(2-/3) 



2n 



r((3/2) 



(2.36) 



(2.37) 



results from the equality r(x + l) = xT(x) for x = (1 - (3)/(2- (3). <P(f3) was 
chosen such that it equals to 1 at the collapse point, namely 



*08) =1 + [C + ln(47r)] ((3-1) + 0((f3 - if). 



(2.38) 



Using the representation (12 .36[) in (|2.35l) . we get the density- fugacity rela- 
tionship which involves the leading hard-core correction: 



n(z, a) 



1-/3 



(2.39) 
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The leading correction oc ^ is marginal in the stability region < (3 < 1. 
It makes the density finite at the collapse point j3 c = 1, 

V&gl = -4 [In £ + C + ln(47r)] . (2.40) 

Note that the density is positive for small £ and diverges for £ — > 0, as it 
should be. The leading correction becomes relevant in the region /? > 1, up 
to p = 3/2 at which <£(/3) diverges. To overcome this divergence, we have to 
account for higher-order terms of the short-distance expansion of the Ursell 
functions in the screening sum rule (|2.34p , in close analogy with [2"3"ll3"Tll3"2"] . 
Such analysis goes beyond the scope of this work. 

It is convenient to pass from the parameter £ to the packing fraction 

■q = no. (2.41) 

Since each particle occupies the length a, it must hold that L/N > cr, i.e. 

2 

rj < 1. Denoting f = n(z, a)/z 2 -f , we have r] — /£. Multiplying both sides 
of Ea. mt bv (1 - 0) f 1 -?. we get 

(1 - p)f-P - 4 ^(Z?)/ 1 ^ - v 1 -?} = 0. (2.42) 



At P = 0, from two solutions of this equation we take the one / = 1 + ^/1 — Aij 
and follow this branch when increasing /3. The density-fugacity relation for 
the packing fractions 10 -3 and 10 -2 are plotted in Figs. 1 and 2 by the 
dashed and dotted curves, respectively. 

For particles with hard cores, in analogy with (|2.14p it holds 

d 

n(z,a)=z—/3P(z,a). (2.43) 
With regard to (|2.35p . we obtain the equation of state 

PP(z,a) =PP(z,0)-2z 2 ^-^, £^0, (2.44) 

where f3P(z, 0) is given by (I2.12[) . This relation can be deduced directly from 
the definition of the grand partition function; the importance of such deriva- 
tion is evident for non-conducting systems which do not exhibit screening. As 

the collapse phenomenon does not involve (++) and ( ) pairs of charges, 

we simplify the model by considering hard cores only between oppositely 
charged particles. In particular, v qq (x, x') — — ln(|x — x'\/ro) for an arbitrary 
distance between particles and 

{hx{\x - x'\/r Q ) if \x - x'\ > a, 
(2.45) 
00 if \x — x'\ < 0". 

With the definition (I2.4j) of the fugacity z, the grand partition function of 
this neutral ID system of size L (which is large enough to neglect boundary 
effects, but not infinite) is given by 

oo 

~ L {z,a) = 1 + ]T z 2N I 2N (a), (2.46) 

N=l 
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where the configuration integrals 

r-L 



hNi<j) = TWJ dxi l dx ' v "L dXN l dx ' N 



N 



iB 2N ({x i },{x' i }) l[ 9{\x i -x' j \-a) 



(2.47) 



involve the Boltzmann weights 



B 2N ({x t },{x' i }) 



Tl(i<:j) = i( x i x j)( x 'i x j) 



Here, 0(x) is the Heaviside theta function: 



9{x) = 



\0i 



if x > 0, 
if x < 0, 



(2.48) 



(2.49) 



which fulfills the obvious relation 9(x) + 0(—x) — 1. Our trick consists in 
substituting each theta function in (I2.47[) by 



9(\xi - x'A — a) = 1 - 0(<r - lac, - x'A). 



(2.50) 



For the lowest-order integral /2(c), we get 

h{o) = I 2 (0)-^— [ 6x x f dx[- 1 - TJ ^9(a - \ Xl - ^ |) 

l i! J Jo Jo Fi *il 

~I 2 (0)-2L/ dxx-' 3 =/ 2 (0)-2L- -. (2.51) 

Jo 1 — P 

For a general integral I 2 n((t), we expand the product of theta functions in 
(ETTT)) as follows 



N 

n 

»>i=i 



iV 



^•i-o-j = 1- E ^-K-^-i) 

AT 

+ E ^-ki-^-l)^-kfc-^|) + 



(2.52) 



The first term on the rhs (unity) gives I 2 n(0)- The second term contains 
N 2 summands, each provides the same contributions because the Boltzmann 
weight (|2.48[) is symmetric with respect to any interchange of two Xi or two 
x\. This is why we can substitute the sum by N 2 9(a — \x\ — x[\). The oppo- 
sitely charged particles at x\ and x[ have to be very close to one another for 
small a, and therefore form an almost neutral entity which decouples from all 
other charges. The integration over x\ and x\ produces the previous factor 
— 2Lcr 1_ ^/(l — while the integration over x-coordinates of all remaining 
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2(N - 1) charges, when multiplied by N 2 /(N \) 2 , implies I 2 (N-i)(0)- The 
summands in the third term on the rhs of (I2.52j) can have i = k or j = /; we 
exclude such terms because they describe the situation of one say (+) charge 
coupled to two (— ) charges and the corresponding power cr 2- ' 3 is marginal in 
comparison with u 1- ^. Then there is [N(N — l)] 2 /2! of equivalent terms of 
type 0(a — \xi — x' x \)9(a — \x 2 — x' 2 \). As before, the couples of particles with 
coordinates (x\,x'i) and (x2,x' 2 ) form neutral entities for small cr, each con- 
tributing by the factor — 2Lcr 1 ~P / (1 — 0). The integration over x-coordinates 
of all remaining 2(N-2) charges, multiplied by [N(N- 1)] 2 /(2!A! 2 ), implies 
7 2 (Ar-2)(0)/2!. Proceeding further in this way we end up with the recurrence 
relation 

I 2N (a) = I 2N (0) - \2Lj=p ) h(N-i) (0) + - { \* L Y=p ) /2 (^- 2 ) (°) + " • " 

(2.53) 

with Iq(0) = 1. Inserting this recurrence into the definition (I2.46|) . we end up 
with 



~ L (z, a) = ~ L (z, 0) exp . (2.54) 

Applying the logarithm to both sides of this equation and using the definition 
of the pressure (|2.8[) . we recover the needed result (|2.44[> . 



3 Thermodynamics of the ordered ID log-gas 

We proceed by the log-gas system with alternating ± charges. To distinguish 
the quantities from the unconstrained case, we add to them the upperscript 
"(ord)", i.e. ordered. As before, we start with pointlike charges in Sect. 13. ll 
and then, in order to pass the collapse point /3 C = 1, we attach to particles a 
small hard core in Sect. 13.21 



3.1 Pointlike particles 

For a neutral configuration of 2N particles constrained to a circle of circum- 
ference L, we consider only ordered sequences with alternating ± charges: 



The grand partition function is expressible as 

oo 

4 OTd )(z) = l+^[(2^/ 2 zi 1 - 



N=l 



2N 



where the configuration integrals 



I. 



(ord) _ 
2N " 



d<gi 
2tt 



V1 M 
2tt 



dy 2 
2tt 



dy 2 
2tt 



(ord) 
L 2N ' 



d<p 



N 



2tt 



xB 2N (Wi},{<p'i}) 



(3.1) 



(3.2) 



N 



2ir 
(3.3) 
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involve the Boltzmann weights (|2 . T[) . We are interested in the thermodynamic 
limit of the pressure p( ord ) and the total particle density n^ ord \ 



pp(° ld \z) = lim -In ~ ( L ° ld) (z), 



i {md \z) = z^-PP^Xz). (3.4) 
oz 



As in the unconstrained case, both ^P(° rd )(z) and n( ord )(z) scale with the 

2 

fugacity like z 2 -f> . 

The ordered configuration integrals (|3.3[) can be compared to their uncon- 
strained counterparts (|2.6[) . The point is that the Boltzmann weights (|2.7p 
are invariant with respect to any interchange of two <fi or two </jJ. Since there 
exist AM possible ways how to order {ifi}, and similarly for {</^}, we can 
choose the special angle constraints 



ipi > tp 2 > ■ ■ ■ > fN-l > <PN, 

and rewrite (12.61) as follows 



ip 1 ><p 2 > ■■■> ip N _ 1 > ip N 



(3.5) 



I-2N — 



d^i 
2tt 



2tt 



d<f2 
2tt 



2tt 



dip 



N 



2tt 



Vjv - 1 d^ 



N 



xB 2N {{i Pi },y l }). 



2tt 
(3.6) 



The alternating configuration space (|3.1[) constitutes a subset of the reduced 
configuration space Q3.5P . Since the Boltzmann weights are positive, we con- 
clude that I 2 n > ^2N or ' equivalently, 5l(z) > S^ ld \z). In the thermo- 
dynamic limit L — > oo, we have the rigorous inequalities 



0P(z) > /3P (ord) (2), n(z) > n (old \z). 



(3.7) 



We take the TBA results from [10] . The expansion parameter in Q3.2p 
is expressible in terms of the length scale Lb in analogy with (|2.9p if the 
substitution z — > rz with r = i/(q— q^ 1 ) = l/[2 sin(7r/3/2)] is made: 



2 sin 



(2^)^ 2 i 1 -^ =r 



r 



i 

2-/3 



(3.8) 



For large L, the grand partition function for the ordered system was found 
to behave as 



~(ord) 



exp 



1^ 



■ tan 



7T/3 



2(2-/3) 



(3.9) 



Eliminating Lb from the two equations and using the formulas (|2.11l) for the 
Gamma functions, the bulk pressure is obtained in the form 



l3p(° rd Xz) 



r 



2-/3 J 



r 



\4-2P J 



m-fi- 



■2 - n 



(3.10) 
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Fig. 3 The ratio of the pressures for the ordered and unconstrained charge systems 
versus the inverse temperature f3. 



For a fixed z, the ratio of this pressure to the pressure for the uncon- 
strained charge system (|2.12[) turns out to be 



p(ord) 
P 



2 sin 



7T/3 



2(2-/3) 



2 sin 



4) 



(3.11) 



This ratio is plotted as a function of the inverse temperature j3 in Fig. 3. It is 
seen that the pressure inequality in (|3.7p is satisfied. Interestingly, the ratio 
takes the minimum value 1 /2 just at the limiting points (3 — > and j3 — > 1. 

In Appendix A, a fusion relation between the grand partition functions of 
the unconstrained and ordered log-gases is used to rederive in an alternative 
way the relationship between the corresponding bulk pressures (|3.11l) . 

Having the bulk pressure (I3.10p . the density- fugacity relationship reads 

as 

2 r (^f) 



,(ord) 



0F(2 -P) r ( 4-3/3 



4-2/3 



r i 



The small-/3 expansion of the rhs of this formula 



.(ord) 



z 2 - 



= I + \{C + I + \n2)P + 0{p 2 ) 



(3.12) 



(3.13) 



is analytic, which is in contrast to the singular result (|2.16[) for the un- 
constrained log-gas. This expansion is checked in Appendix B, using pair 
distributions of the free (/3 = 0) log-gas with alternating ± charges pQ. For 
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a fixed z and in the limit the term T((l - /3)/(2 - 0)) - 1/(1 - (3) 

implies the collapse singularity 

« (ord) ~ T^T (3-14) 

The adequacy of this term, in the sense that it compensates exactly the singu- 
larity of the hard-core contribution at the collapse point, will be documented 
in the next part. The plot of the density- fugacity relationship for ordered 
pointlike particles (|3.12p . valid up to /3 C = 1, is represented in Fig. 4 by the 
solid curve. The curve is monotonously increasing in the whole /3-interval 
[0, 1], which is in contrast to Fig. 2 for the unconstrained system. 




Fig. 4 The density-fugacity relationship versus the inverse temperature /3 for the 
ordered log-gas. The solid curve describes pointlike charges with the packing frac- 
tion r) = 0. The dashed and dotted curves correspond to -q = 10 -3 and 10 -2 , 
respectively. 



The virial equation of state takes the usual form 
/?P (ord) = (l-^j n^ d \ 

The (dimensionless) Helmholtz free energy per particle reads 



f^ifl, n) = ln(An) - £ ln(r„n) + £ ( 1 - £ ) In tt 



P 



l ~2 



l-ln(l-| 



lnT 



1- /3 

2- /3 



lnT 



_3/T 
~20, 



(3.15) 



(3.16) 
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For the excess specific heat, we obtain 



c (ord)cx p2 



k B 2(2 - p) 



V \2-p) V \A-2p))\ 



02 + ^'(i-ri. (3.7, 



2(2-/3) 4 r V 2 
Close to the collapse point /3 — > 1 — , it exhibits the singular behavior 
ci ord)cx 1 3 



fc B 2(1 -fl a 2(1-/3) f °W< ^ 
which is exactly the same as in the unconstrained case, see Eq. (|2.28p . 

3.2 Particles with hard cores 

We attach to particles hard cores of diameter a and derive the equation of 
state directly from the definition of the grand partition function, in close 
analogy with the procedure for the unconstrained system starting from Eq. 
(I2~il) . 

For the ordered charge system of large size L, the grand partition function 
is defined by 



S<^W) = 1 + £ ^4r>)> (3-19) 

N=l 

where the configuration integrals 

(Old) f L f X N-l f x N 

^2N ( a ) = / / dx[ / dx 2 J dx' 2 ■ ■ ■ I dxN / da;^ 

Jo Jo Jo Jo Jo Jo 

N N-l 

xB 2N ({x i },{x' i })l[6(x i -x' i -a) ]J 6(x\ - x l+1 - a) (3.20) 

i=l i=l 

involve the Boltzmann weights (|2.48l) . As before, we substitute each Heaviside 
theta function in (|3.20p via (|2.50p . For the lowest-order integral, we get 

4° r <V)=4° rd) (0)- f L dxi f 1 dx' 1] 1 -—0(a-x 1 +x[) 

Jo Jo Fi — x ir 

~/ 2 (0)-L / dxx-P = I 2 (0)-L^- (3.21) 
Jo 1 P 



For a general integral I 2 °^ in p.20p . we expand the product of theta func- 
tions written as in (|2.50p . The first term (unity) implies I 2 °^ d \o). The second 
term contains (2N — 1) theta functions; each theta function forces two op- 
positely charged nearest-neighbor particles to be very close to one another 
for small a, and therefore to form an almost neutral entity which decouples 
from all other charges. The integration over distance between these nearest 
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neighbors produces the factor — er 1_,9 /(l — /3). The integration over the re- 
coordinate of the neutral entity, which exhibits all possible orderings with 
respect to a given configuration of other "active" charges, implies L. The re- 
maining 2(N — 1) charges contribute by ^(jv— so ^ ne total contribution 
is [— Ler 1 ~' 9 /(1 — /3)]^2(jv-i) (^)- I n the case of k neutral entities, there is an 
indistinguishability factor 1/fc! attacl 
end up with the recurrence relation 



indistinguishability factor 1/fc! attached to [—La 1 ^ 13 /(l~/3)} k I^l k ^(0). We 



4°; d V) =4°; d) (o)-(^) 4(^(0)+^ (if^) 2 4(; d i 2) (o)+--- 

(3.22) 

with Ig° ld \o) = 1. Inserting this recurrence into the definition (|3.19p . we end 
up with 

~f d) (z,a) = ~f d >(z,0)exp (-* 2 f^) • (3-23) 

Applying the logarithm to both sides of this equation and in the limit L — > oo, 
we get the bulk pressure of the ordered system of charges with hard cores: 

/3F(° rd )(z,a) = /3P( OTd )(z,0)-z 2 ^, £->0, (3.24) 

where /3P^(z,0) is given by p. 101) . The corresponding density- fugacity 
relationship reads as 

n^ d \ Zl cj) = n (md \z,0)-2z 2 — -, £ ->■ 0. (3.25) 

1 p 

Although these results were derived in the stability region < (3 < 1, we 
assume that they can be analytically continued beyond the collapse point 
/3 C = 1- 

To provide a convenient representation of p. 251) , we express the density- 
fugacity relation for pointlike particles (13. 12)) as follows 

n (ord) (z, 0) = _ o z ^ ' (3 ' 26) 
i p 



where the function 



<Z>( OTd )(/3) = 



1 F \ — 



3-2/3 



r„-§ 



(3.27) 



was chosen such that it equals to 1 at the collapse point, namely 

^(ord)^ = i + (C + In tt) (/3 - 1) + O(03 - l) 2 ). (3.28) 
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Using the representation (|3.26p in (|3.25l) . we get the density- fugacity rela- 
tionship which involves the leading hard-core correction: 



,(ord) 



0,ct) 



(3.29) 



The leading correction oc P makes the density finite at the collapse point 

/9c = 1, 

p?± = -2 (In £ + C + In tt) . (3.30) 

The formalism holds up to f3 = 3/2 at which & ovA Xfi) diverges. 

2 

It is convenient to pass from the parameter £ = az 2 -f to the packing 



tract 



ion r\ 



= r i (° rd ) 



a. Denoting g 



= r).( ord )f*. nM?.— 



plying both sides of Eq. (pT2"9"]) by (1 - /3).g 1_/3 , we get 

(1 _ p)g*-P _ 2 U OTd )(/3) 3 1 -' 3 - if-? 



(z, <j)/z 2 -p , we have rj = g£. Multi- 



0. 



(3.31) 



At P = 0, we take the solution g = (1 + \/l — 877) /2 and follow this branch 
with increasing /?. The density-fugacity relation for the packing fractions 10~ 3 
and 10 -2 are plotted in Fig. 4 by the dashed and dotted curves, respectively. 
The functions are monotonously increasing in the whole interval of j3, which 
is in contrast to the unconstrained system (see Fig. 2). 



4 Conclusion 

The aim of this work was to derive and analyze the exact thermodynamics 
of the two-component log-gas formulated on an infinite line with alternating 
± charges. 

The thermodynamics was first solved for the case of pointlike particles, up 
to the collapse point /3 C = 1, by using the TBA results. The obtained density- 
fugacity relationship (|3.12l) turns out to be analytic in the high-temperature 
region, see the small-/? expansion (I3.13[) and its check in the Appendix B. 
Moreover, as is shown in Fig. 4 by the solid curve, the density-fugacity ratio 
n (oid) j z 2/{2-p) j g a mono tonously increasing function of the inverse temper- 
ature (3. This is in contrast to the unconstrained log-gas with a non-analytic 
small-/? expansion and a non-monotonous plot of n/z 2 '( 2 ~^ versus /?, see the 
solid curve in Fig. 2. For a fixed fugacity z, the validity of the inequality (13. 7|) 
for the pressures of the unconstrained and ordered log-gases was confirmed 
by finding the exact ratio p( ord )/P (Fig. 3); the minimum 1/2 is reached just 
at the limiting points /? — > and j3 — > 1. 

The inclusion of hard cores around the particles was a more complicated 
task than in the unconstrained case since we could not apply the screening 
sum rule for pair correlation functions valid for a conductor. The derivation 
of the leading hard-core correction to the pressure, outlined in Sect. 3.2, 
was based on the explicit definition of the grand partition function and a 
substitution of Heaviside theta functions due to the presence of hard cores. 
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The final result for the bulk pressure (|3 . 24[) and the corresponding density- 
fugacity relationship (I3.25[) pass the collapse test, in the sense that the hard- 
core correction compensates exactly the singularity of the pointlike model 
and leads to a finite pressure and particle density at p = 1. 

Since the log-gas is by the structure much simpler that the 2D Coulomb 
gas, it might be possible to get all relevant hard-core corrections which com- 
pensate an infinite series of singularities of the pointlike result and allow us 
to go up to /3kt = 2, in the spirit of Refs. [231131 1132] . This problem is left for 
future. 



Appendix A 



Here, we rederive in an alternative way the relationship between the bulk 
pressures of the ordered and unconstrained log-gases (|3.11l) . 

We use the fusion relation between the grand partition functions of the 
unconstrained and ordered log-gases [T2TIT5] : 



iold \( q - q - 1 )z) = 



S L (qz) + ^. L (q 1 z) 



q = c^/ 2 . 



(A.l) 



2S L (z) 

We know that for large L the grand partition functions behave as 

pP {md \z)L\ , (3P {md \z) = a {md \p)z^e (A.2) 



-(ord) 



and 



(z) ~ exp 



S L {z) ~ exp \PP{z)L] , 0P(z) = a{P)z-*=v . 



(A.3) 



The lhs of Eq. flgjj reads 



exp|-La(° rd )(/3) 

La (ord) (/3) 



2z sin 



7T/3 



X COS ■ 



2z sin — — 



sm 



2- J 



7T/3 



2(2-/3) 
nf3 



2(2-/3) 



(A.4) 



The rhs of Eq. (|A.1[) can be expressed as 

7T/3 



exp < La{j3)z l 



2-(3 



Lai^z 2 -/ 3 sin 



irf3 
2^3 



(A.5) 

The equality of separately exponential and cosine terms in Eqs. (|A.4[) and 

7T/3 



(|A.5|) leads to the only relation 

(ord)^ ^ p(ord) 



a(P) P 
in agreement with (|3.11l) 



= 2 sin 



2(2-/3) 



2 sin 



(A.6) 
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Appendix B 

Here, we construct the small-/3 expansion of the pressure and the particle 
density for the ordered log-gas. 

For (3 = 0, the Boltzmann factor E>2N ({<Pi} > Wi}) = 1- The ordered 
integrals Q3.3P are simply given by J^' d) = 1/(2 AT)! as there are just (2N)\ 
ways how to order the angle variables and each ordering implies the same 
contribution. This leads to 



2N 

1+ £w-°> p(zi) - (R1) 

N=l y ' 



Consequently, in the lowest order, 

f3Pt d \z) = z, n^ d) (z)=z. (B.2) 

For a given configuration of particles {j} with charges {qj} at positions 
{xj}, the interaction energy is given by 

E = - ^2qjq k (-ln\xj - x k \) 



2 



\ f\x f dx'(-hi\x-x'\) 

1 J ° J ° q,q'=±l 



qq'h qq >(x,x'), (B.3) 



where the microscopic quantity 



h gg r(x,x') = ),6 q!q .5(x - Xj)S q ' <gh S(x' - x k ). (B.4) 

Expanding the Boltzmann factor in (3, e~^ E ~ 1- (3E, and using the cumu- 
lant expansion, we get the leading /3-correction to the pressure, 

f$p^\ z )=f} P ^\ z )-P^, (B.5) 

where the thermal averaging (■ ■ -)q is over the system of non-interacting (f3 — 
0) ordered charges. For an infinite system L — > cxd, the two-body distributions 
n qq >(x,x') = n qq i(\x — x'\) = (h qq > {x, x'))q have been calculated in [T|: 



l qq\ 



,(*) = (f) 2 (l " e" 2 " 1 * 1 ) . n q ,. q {x) = [I)' (l + e-N) . (B , 
Considering the lowest order for n from (|B.2[) . we obtain 



j3P[ md \z) = z + — j dx \n\x\[n ++ (x) + n (x) — (x) — n \-{x)] 

^ J -00 

= z + ^z(C + \n2) + ^z\nz. (B.7) 
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The corresponding density of particles 

n[ md \z) = z+^z(C+l + ln2) + ^lnz (B.8) 

is consistent with the small-/? expansion (13.131) . where z^~$ ~ z[l+(/3/2) In z\. 

To find the next order of the small-/3 expansion we need to know four- 
body distributions for the system of non-interacting ordered charges, and so 
on. 
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